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The role of intrinsic fex in inclusive and semi-inclusive Deep Inelastic Scattering processes {£p ^ 
IhX) is studied with exact kinematics within QCD parton model at leading order; the dependence of 
the unpolarized cross section on the azimuthal angle between the leptonic and the hadron production 
planes (Cahn effect) is compared with data and used to estimate the average values of fc^ both in 
quark distribution and fragmentation functions. The resulting picture is applied to the description of 
the weighted single spin asymmetry ^^y'^""'*^' recently measured by the HERMES collaboration 
at DESY; this allows to extract some simple models for the quark Sivers functions. These are 
compared with the Sivers functions which succeed in describing the data on transverse single spin 
asymmetries in p^p — » ttX processes; the two sets of functions are not inconsistent. The extracted 
Sivers functions give predictions for the COMPASS measurement of Al^^'^^'^'^^^ in agreement with 
recent preliminary data, while their contribution to HERMES AJ)"**" is computed and found to be 
small. Predictions for ylj}^''*'^"'*^^' for kaon production at HERMES are also given. 

PACS numbers: 13.88.+C, 13.60.-r, 13.15.+g, 13.85.Ni 



I. INTRODUCTION 



It is becoming increasingly clear that unintegrated parton distribution and fragmentation functions play a significant 
role in many physical processes and that they are more fundamental objects than the usual integrated versions; 
intrinsic fcj^ originates both from partonic confinement and from basic QCD evolution P| and often cannot be ignored 
in perturbative QCD hard processes and in soft non perturbative physics. For example, this was already pointed out 
in Refs. concerning the computation of the unpolarized cross section for inclusive hard scattering processes like 
pp h X at intermediate energy values. Similarly, in semi-inclusive Deep Inelastic Scattering (SIDIS) processes, 
ip^ £hX, the intrinsic partonic motion results in an azimuthal dependence of the produced hadron A 
full consistent treatment of several inclusive pp processes with all intrinsic motions, in different kinematical regions, 
has been recently discussed in Ref. 0- 

While the role of intrinsic k± can be important in unpolarized processes, it becomes crucial for the explanation of 
many single spin effects recently observed and still under active investigation in several ongoing experiments; spin and 
kj_ dependences can couple in parton distributions and fragmentations Q , thus giving origin to unexpected effects in 
polarization observables. One such example is the azimuthal asymmetry observed by the HERMES collaboration in 
the scattering of unpolarized leptons off polarized protons 0, ^3 • Another striking case is the observation of large 
transverse single spin asymmetries (SSA) in p — > tt X processes 0, . 

We consider here in a consistent approach the role of parton intrinsic motion in inclusive and semi-inclusive DIS 
processes within QCD parton model at leading order; the kinematics of intrinsic fc_L is fully taken into account in 
quark distribution functions, in the elementary processes and in the quark fragmentation process. This induces several 
corrections of order k±/Q or higher, which are exactly computed: however, we do not consider similar corrections 
which might originate from higher- twist distribution and fragmentation functions [l^ . The average values of k±_ for 
quarks inside protons, and for final hadrons inside the fragmenting quark jet, are fixed by a comparison with data on 
the dependence of the unpolarized cross section on the azimuthal angle between the leptonic and the hadronic planes. 

Such values are then used to compute the SSA for ^ in X processes, which would be zero without any intrinsic 
motion. We concentrate on the Sivers mechanism that is the spin and k± dependence in the distribution of 

unpolarized quarks inside a transversel y po larized proton; it can be isolated by stud ying the weighted SSA A^'^'^" 
|15j . recently measured by HERMES [13 and, still preliminarily, by COMPASS [ia collaborations. Although the 
data are still scarce, with large errors, a first qualitative estimate of the quark Sivers functions can be obtained; the 
information gathered from HERMES data is in agreement with the preliminary COMPASS data. The contribution 
of these functions to the weighted longitudinal SSA A^"^"^" is computed and shown to be negligible: the measured 



FIG. 1: Kinematics of DIS processes. 



-^iiif^ can be equally originated not only by the Sivers mechanism, but also by the Collins effect occurring 
in the fragmentation of a transversely polarized quark, and by higher-twist contributions. 

A similar analysis of SSA in p — tt ^processes, with a separate study of the Sivers and the Collins contributions, 
has been performed respectively in Rcfs. Q and |l8j . with the conclusion that the Sivers mechanism alone can explain 
the data [llj, while the Collins mechanism is strongly suppressed. Explicit expressions of the quark Sivers functions, 
as obtained from data on SSA in p n X processes, are given in ref. 0. They are qualitatively similar to the 
Sivers functions obtained here from SIDIS processes. Let us notice that the universality of the Sivers function is an 
important open issue; it has been proven that the quark Sivers functions in SIDIS and in Drell-Yan processes must be 
opposite |l9| . while no definite conclusion has been theoretically reached concerning the relation between the Sivers 
function in SIDIS and pp processes |23| • 



II. DEFINITIONS AND KINEMATICS 



Let us start from the kinematics of Deep Inelastic Scattering processes in the j*p cm. frame, as shown in Fig. ^ 
We take the photon and the proton colliding along the z axis with momenta qf and P respectively; the leptonic plane 
coincides with the x-z plane (following the so-called "Trento conventions" |21|). We adopt the usual DIS variables 
(neglecting the lepton mass): 

2P-q W^ + Q^-ml ^ P-l xu{s-ml)' ^' 

If one neglects also the proton mass nip the four-momenta involved can be written as: 

i = E{l,sm6,0,cose) {E,£) 

e = l-q 

P = Po(l, 0,0,-1) , 
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where 



cos 9 



Q^ + E{W-QyW) l + (y-2): 



E{W + Q^/W) 



E ■ 



1 



2W 



(2) 



At leading QCD order the lepton scatters off a quark and, taking intrinsic motion into account, the initial and final 
quark four- moment a are given by: 



xPo 



-xPo 



ixPoJ 



k' = k + q 



(3) 



where x ~ k^ / is the light-cone fraction of the proton momentum carried by the parton and k± = k± (cos ip, sin ip, 0) 
is the parton transverse momentum, with k± = \k±\. 

The elementary Mandelstam variables s = {I + fc)^, i = {I — l')'^, and u ^ {k — I'Y are then given by: 



xs-21- kj_ 



(4) 



-x{ s-—] + 2£-ki 



xQ2 



The on-shell condition for the final quark 



= 2(7 • fc - 0^ = s + t + u = 



(5) 



implies 



(6) 



Notice that when terms 0{k\/Q'^) are neglected in the above equations one recovers the usual relations x = Xb 
and fe = XbP + k^] however, a significant dependence on the azimuthal angle remains, at 0{k±/Q), in the partonic 
Mandelstam invariants s and u, via 



£ ■ k± = Ek± sin 9 cos ip . 



(7) 



One explicitly has 



y 
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1 - 4— y^l-y cos tp 



O 



kl 



^ fc_L cos if 



o 



(8) 
(9) 



A. DIS cross section 



Let us recall the expression of the cross section for Deep Inelastic Scattering processes, £p ^ £ X, in the framework 
of the QCD parton model with the inclusion of intrinsic k± One starts from 

with 

Lfj,^ = 2{lfj.ll + - gf^^l ■ I') , (11) 
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and 

^""^ = E/ dxd^k^ f,{x,k^)w^''' , (12) 

where fq{x, kj_) is the number density of quarks of flavor q inside the initial hadron. carrying a transverse momentum 
fcj^ and a hght-cone fraction x of the proton momentum. The elementary quark tensor is given by 



w 



2el ( + k'^k. ~9^.^) 5{2q ■ k ~ Q^) , (13) 



flu 
so that 

L^, w''" = 2 el S{2q ■ k ~ Q^) (s^ + u^) . (14) 
In the coUinear case, fq{x,k±) — fq{x) 6'^{k±), Eqs. ((Tm) - ((Tl|l simply result into 

dx, dQ^ - 2^ fii^B)^^ , I15j 

where da''^^'* is the cross section for the elementary lepton-quark scattering, 



9 :j;2s2 Q 

In the general non coUinear case one obtains instead 
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where x is given in Eq. lO and the function J is given by 



E / Mx,k^) J(a;„Q2^fci) , (17) 



Notice that J = 1 in the coUinear case; the elementary da in Eq. H17() is the same as in Eq. 1)16(1 with the Mandelstam 
variables of Eqs. which depend on k±. 

If one could detect the final quark - for example by reconstructing the current fragmentation jet ~ the DIS cross 
section could be written as 

and one could test the azimuthal dependence of the cross section on the angle ^) between the leptonic and the jet plane, 
Fig. 1. Such a dependence, resulting from + , Eqs. © and ®, was suggested by Cahn Q, in semi-inclusive 
DIS, assuming that the fragmentation process q — > hX is essentially coUinear and the direction of the final detected 
hadron is close to that of the quark. Smearing effects due to the fragmentation process were also taken into account 
[M^Ull' In the next subsection we shall consider SIDIS processes, fully taking into account the intrinsic motion and 
the angular dependence in the fragmentation process; we shall see that indeed a dependence on an azimuthal angle 
remains, which allows to obtain an estimate of the average transverse momentum in distribution and fragmentation 
functions. 



B. SIDIS cross section 



Let us consider semi-inclusive DIS processes, — > thX^ in the 7*p cm. frame and in the kinematic regime 
in which Pt ~ Aqcd — where Pt = \Pt\ is the final hadron transverse momentum. In this region leading 
order elementary processes, Iq^ iq.^ are dominating: the soft Pt of the detected hadron is mainly originating from 
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FIG. 2: Kinematics of the fragmentation process. 



quark intrinsic motion |5LI3.|24|, rather than from higher order pQCD interactions, which, instead, would dominantly 
produce large Pt hadrons |2,'ll2(l . 

We adopt as our starting point a factorized scheme, introducing - within the leading-order parton model with 
leading-twist distribution and fragmentation functions - the exact k± kinematics; this induces corrections of the 
0{k±/Q) or higher, which we wish to explore and compare with experiments. A formal QCD factorization, in the 
same low transverse momentum region and at leading order in kj_/Q, has been recently proved p^ . 

In the factorization scheme, assuming an independent fragmentation process, the SIDIS cross section for the pro- 
duction of a hadron h inside the jet originated from a final quark with transverse momentum fej^ can be written 
as 

j^7^lp^e+jet+h+X 



dxB dQ^ d'^k± dz d'^p^ 



fq{x,kj_) 



da^ 
~dQ^ 



(20) 



where D^{z^ p^_) is the number density of hadrons h resulting from the fragmentation of the final parton q, normalized 
so that 

Jdzd^p^D^^iz,p^)^(Nh), (21) 

where (Nh) is the average multiplicity of hadron h in the current fragmentation region of quark q and 

jd^p^D'^iz,p^) = D'^iz). (22) 

p^ is the transverse momentum of the hadron h with respect to the direction k' of the fragmenting quark and 
z = Pf^/k^ is the light-cone fraction of the quark momentum carried by the resulting hadron in the (a;, y, z)-system 
(see Figs. |21and|31). These natural variables for the fragmentation process can be written as: 



k 



2k'" 

p^ = Ph-{Ph-k)k , 
where P,, = (P^, Pt, Pf,) and k' = k + q = (fc'°, k^,k'^), with: 



(23) 
(24) 



^,0 _ ^x- 2x, 



k- = ^ 



Ife'l^ = kj 



1 3/ 



l~ Xb 

{k'^f 



(25) 



The above equations allow us to express, for each value of Xg and Q , the fragmentation variables z and p^ in 



terms of the usual observed hadronic variables Pt and Zh ^ {P ■ Ph)/{P ■ q) ~ {P^ + Ph)/^- O^*^ finds 



1 



2fc'° 



ZhW 



p2 



2 2z,,I^ 
k^Pr 2x, 



fe, 



Ifc'l 



k^ 



ZhW 



2zhW 



1 



(26) 
(27) 
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FIG. 3: Three dimensional kinematics of the SIDIS process. 




and 



Pr-k^ + Pjk'^ 3 Pr-k^+Pjk'^ ,3.. 
- k^,P, k j (28) 

-Zhk^ + o(^) (29) 



where fc'O, and \k'\ are given in Eqs. ^ and P^ = {zh W)/2 - P^/{2zh W). 

Eqs. (|26|l and H28|) allow us to describe the fragmentation process in terms of the variables (zh, Pt)- 

dz (fp^ = dzh d^Pr — , (30) 

Zh 

so that, finally, the SIDIS cross section (|20|l can be written in terms of physical observables as: 



This is an exact expression at all orders in {k±/Q); x is given in Eq. © and the full expressions of z and in 
terms of XB,Q'^,k±,Zh and Pt can be derived from Eqs. (|25|l . H26() and (|28|l . Notice that, in the physical variables 
a;^ and z^i, the x — z factorization of Eq. H2U|) is lost, even in our simple parton model treatment; it can be recovered 
at 0{kjQ) (see Eq. ^ below). 

Let us now consider again the issue discussed at the end of Section fll Al concerning the azimuthal dependence of the 
cross section, by comparing Eqs. H19() and H31() . The former equation describes the cross section for jet production and 
depends, as we explained, on the azimuthal angle (f, that is on the azimuthal angle of the intrinsic fcj^ of the quark in 
the proton. Such a dependence is integrated over in Eq. which describes the cross section for the production of a 
hadron, resulting from the non coUinear fragmentation of the quark. Therefore, there cannot be any tp dependence in 
this cross section. However, due to relations and the integration over fc_L at fixed Pt = /^(cos sin(/)/i, 0) 
introduces a dependence on the azimuthal angle of the produced hadron h, that is the angle between the leptonic 
and the hadronic plane, Fig. |31 This azimuthal dependence remains in the SIDIS cross section and will be studied in 
the next Section (see also Appendix A). 
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It is instructive, and often quite accurate, to consider the above equations in the much simpler hmit in which only 
terms of 0{k±/Q) are retained. In such and one obtains: 



dxs dQ^ dzh d'^Pr 
where ~ Pt — Zh k±, Eq. (|29|l . and 



D'^{zh,PA.), 



(32) 



(^1 + (1 - yf - 4:^(2 - y)v/T^cos^^ 



(33) 



In what follows we assume, both for the parton densities and the fragmentation functions, the usual factorization 
between the intrinsic transverse momentum and the light-cone fraction dependences, with a Gaussian k± dependence, 
that is: 



fq{x,kj_) = fg{x) 



1 



n{kl) 



and 



so that 



and 



D'iiz,p^) ^ D'liz) 



-Pl/ipl) 



d?k±_ fq{x,kj_) = fg{x) 



I d^p^ D'^iz,p^)^D'^q{z). 



(34) 



(35) 



(36) 



(37) 



With the above expressions of fq{x, kj_) and D'^{z,p±) the d^k± integration in Eq. 1)32(1 can be performed analyti- 
cally, with the result, valid up to 0{k±/Q): 



dxB dQ^ dzh d^Pr 



where 



, {2-y)^T^{kl)zhPT . 

mQ ^\ 



{P^) = {pi) + zi{ki) 



1 



i + (i-y)' 

PtKPt) 



(38) 



(39) 



This approximate result illustrates very clearly the origin of the dependence of the unpolarized SIDIS cross section 
on the azimuthal angle 0^ . As observed first by Cahn |3. l23l| . such a dependence is related to the parton intrinsic 
motion and it vanishes when fcj^ = 0. Having also taken into account the intrinsic motion in the fragmentation process, 
Eq. H38() also depends on (p^), via the quantity (P^) defined in Eq. (|39|l . 

As we said, the above results hold in the small P^ — Aqcd — ^± region, where corrections ©(fc^/Q^) are expected 
to be small. As we shall see in the next Section the numerical results obtained from Eq. (|31|l or from Eq. H38|l are 
indeed very close. 



III. CAHN EFFECT IN UNPOLARIZED SIDIS 



We wish to obtain experimental information on the average intrinsic motions. Our strategy is that of trying to 
describe several sets of experimental data, which explicitly measure the dependence of the SIDIS unpolarized cross 
section on the azimuthal angle (ph between the lepton plane and the hadron production plane, and on the transverse 



8 




FIG. 4: Fits to the cos (ph dependence of the cross section: the dashed line is obtained with exact kinematics, while the solid 
bold line includes terms up to 0{k±/Q) only. The data are from Ref |^. 



momentum of the detected hadron Px (see Fig. 3); we do that by exploiting Eq. H31|) or H38|) and by fixing the values 
of (fc^) and which best describe the data. 

The (j>h dependence of the SfDIS cross section, for the production of charged hadrons, has been extensively studied 
by the EMC collaboration in the scattering of 280 GcV muons against a hydrogen target HQ. The shape of the 
differential cross section 

dx^ dQ^ dzh Pt dPr 3 — ^ (40) 



d(f)h J dxs dQ^ dzh d'^Px 

is studied as a function of (ph- 

The integration covers the x^, Q^, Zh and Pt regions consistent with the experimental cuts [^: 

xf > 0.1 Pt > 0.2 GcV/c y < 0.8 > A {GcV/cf , (41) 

where xf ~ 2Pl/W and Pl is the longitudinal momentum of the produced hadron relative to the virtual photon, 
according to Fig. 3. 

In our analysis we adopt the MRST 2001 (LO) (2^ parton density functions and the fragmentation functions into 
charged hadrons from Ref. ■ 

Fig. ^ shows our fits to two sets of data 0|, corresponding to two different ranges of xf- The solid bold hue 
shows the result we obtain by taking into account only 0{k±/Q) contributions, Eq. (|38|l . whereas the dashed line 
corresponds to the exact result at all orders in k±/Q, Eq. 131|l . The (— cos^;,) behaviour, explicit in Eq. (|38|l . is 
indeed shown by the data. A possible positive contribution from a cos(20ft) term seems to be visible at small values 
of (pfi (dashed line). 

Other interesting EMC data concern the P^ dependence of the cross section for fip and fid scattering at 
incident beam energy between 100 and 280 GeV; also these data are strictly related to the values of (A:^) and (p^). 
The quantity measured is given by 

I da 1 /■ , ,^2 , d^a'^P'-ehx 

dcpudxedQ dzh - — — „ , (42) 



(JDis dP^ 2 aois J dxg dQ'^ dzh d'^Px ' 

where aujs is the integrated DIS cross section from Eq. H10() . In the integration of Eq. H42() the following experimental 
cuts have been imposed (see Ref. for further details): 

Q2 > 5 (GeV/c)2 , ^ 90 GeV^ , > 5 GeV 

0.1<z,i<0.9, 0.2<y<0.8. (43) 

Fig.Elshows the comparison of our results with EMC data [sjl for different ranges of z^. The solid and dashed lines, 
which are here basically indistinguishable, are the results of our fits at first order and at all orders in k±/Q respectively. 
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0.1 < Zf, < 0.2 



0.2 < Zf, < 0.4 
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0.4<z^< 1.0 
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FIG. 5: The normalized cross section da/dP^: the dashed hne is obtained with exact kinematics, while the solid bold line 
includes terms up to 0{k±/Q) only. The shadowed region corresponds to varying the parameters {k']_) and (pi)) given in Eq. 
lH7|l of the text, by 20%. The data are from Ref. [s^- 



The shadowed region is spanned by varying the parameters (/c^) and (p^) by 20% and shows the sensitivity of our 
results on these parameters. The figure clearly show that, as expected, our LO approach is vahd for Pt values up to 
about 1 GeV/c. At higher values NLO contributions from "f*q gq and "f*g qq processes have to be taken into 
account. 

Useful data on the iph and Pt dependence were also found by the FNAL E665 collaboration 0| in np and fid 
interactions at 490 GeV. The quantity studied is 



(cos 



/ dxB dQ"^ dzh d^Pr cos (j)h d^a 



J dxs dQ-^ dz/i d'^Pr d^a 



where d^a denotes the fully differential cross section 



d^a 



dxg dQ^ dzh d^Pr ' 



(44) 



(45) 



and where the integral over Pt runs from 
integration region in Eq. H44|) is defined by: 



cut pmax 



2.5 GeV/c. According to the experimental setup |2J|, the 



Q2 > 3 (GcV/c)2 , 300 <W^ < 900 GeV^ , 
60 < 1/ < 500 GeV , > 8 GeV , 0.1 <j/< 0.85. 



(46) 



Fig. El shows the data |2J| as a function of P™* compared with the results of our calculations; again, the solid bold 
line corresponds to the result we obtain by taking into account only 0{k±/Q) terms, Eq. (|38() . whereas the dashed line 
corresponds to the exact kinematics, Eq. (|31|l . We also show the EMC data on the xp dependence of {cos(f>h) /wi{y) 
; they compare well with our results obtained by using Eq. (in the right panel of Fig. the theoretical curve 
corresponds to the calculation of (cos0/i)/ {'Wi{y))) with the experimental cuts of Ref. 0| and no integration over Zh 
which is expressed in terms of xp- The shadowed region is obtained by varying the parameters (fc^) and (p^) by 
20%. Once more, Fi g. El clearly shows that our calculation is valid for Pt values up to about 1 GeV/c, where NLO 
pQCD contributions [251 l26l| must be taken into account and our simple LO treatment can no longer be applied. 

All sets of data described above depend crucially on the intrinsic motion in distribution and fragmentation functions; 
their combined analysis leads to the following best values of the parameters: 



(fci) = 0.25 (GeV/c)2 



(pi) = 0.20 (GeV/c)^ 



(47) 



One should notice that the above values have been derived from sets of data collected at different energy, Xb, 
and z^ ranges, looking at the combined production of all charged hadrons in SIDIS processes, and assuming 
constant values of (fci) and (pi), which allow analytical integration up to 0{k±/Q); these values are also assumed 
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FIG. 6: {cos(t>h), as given by Eq. (1441 . as a function of P^"*, left panel, and xf, right panel: the dashed line is obtained with 
exact kinematics, while the solid bold line includes terms up to 0{k±/Q) only. The shadowed region corresponds to varying 
the parameters {k']_) and (px)i given in Eq. 14711 of the text, by 20%. wi{y) = ^^^+(1^^^ ■ data are from Ref. |2^. left 
panel, and Ref. |^, right panel. 



to be independent of the (light) quark flavor. A more refined analysis, introducing for example x and z dependences, 
would require the introduction of new unknown functions. At this stage, we stick to the constant best values (|47f) 
which, together with Eqs. H34|) and H35|) . can only be considered as a consistent simple estimate and a convenient 
parametrization of the true intrinsic motion of quarks in nucleons and of hadrons in jets, supported by the available 
experimental information. In the next Section we adopt such a picture for the computation of pion and kaon SSA in 
SIDIS polarized processes, ^ in X and £p'^ ^ £K X. 



IV. SIVERS EFFECT IN POLARIZED SIDIS 

In this Section we consider the single spin asymmetry Aut measured in ep^ enX processes by the HERMES 
collaboration at DESY, using a transversely polarized proton target [lol |. Our aim is that of obtaining information 
on the quark Sivers functions, which can be isolated and directly accessed by studying the weighted transverse spin 
asymmetry A^!^^*^" "^^^ ^M- These functions arc then compared with those obtained from the study of SSA in p^p 
processes. They are also used to estimate the contribution of the Sivers mechanism alone to the weighted SSA 
^[/l"^" ' measured by the HERMES collaboration in the lepton scattering off a longitudinally polarized proton target 
0. Moreover, expectations for A^!^^*^" in the COMPASS kinematical regions are computed and compared with 

preliminary COMPASS data [Hi and predictions for kaon asymmetries A'J^^'^'^"'^^^ at HERMES are given. 

Let us recall the origin of the SSA in SIDIS, as originated by the Sivers mechanism fl^ . The unpolarized quark (and 
gluon) distributions inside a transversely polarized proton (generically denoted by , with p^ denoting the opposite 
polarization state) can be written as: 

[x, k^) = fg/p{x, k±) + ^ ^^fq/pr {x, k^) S^-iPx fex) , (48) 

where P and are respectively the proton momentum and transverse polarization vector, and k± is the parton 
transverse momentum; transverse refers to the proton direction. Eq. H48|l implies 

fq/pr{x,kj_) + fq/pi{x,k±) = 2fq/p{x,k±) , 

/9/pT (X, fc_L) - fq/pl (X, fcj^) = ^'^fq/pl (x, kj_) -{P X fe_L) , (49) 

where fq/p{x, k±) is the unpolarized parton density and fq/pt (x, k±) is referred to as the Sivers function. Notice 

that, as requested by parity invariance, the scalar quantity ■ {P x kj_) singles out the polarization component 
perpendicular to the P — k± plane. For a proton moving along — z and a generic transverse polarization vector 



11 



= \S^ \ (cos (j)s , sin (f>s , 0) (sec Fig. I^J one has: 

S^-{Px fcj.) = \S^\ sin((^ - (j)s) = l^^l sin0s^. , (50) 

where [ip — (f>s) = ipsiv is the Sivers angle. 

The cross section for the scattering of an unpolarizcd lepton off a polarized proton, in the configuration of Fig. O 
can then simply be written as, see Eq. (|31(l . 



where (f'a^ stands for 

dxB dQ"^ dzhCPPr d(j)s ' ^ 

and where the (jjg dependence is contained in fq^pt{x,k±). The (j)s dependence originates from the fact that, with 
transversely polarized protons, the cross section depends also on the angle between the polarization vector and the 
leptonic plane; in the configuration of Fig. 3 this is simply (jjs — (j^f — 'Psi having chosen (p^i to be zero, with <ps 
varying event by event. In actual experiments variables are measured in a different frame (for example the laboratory 
frame where the lepton moves along the z-axis and the proton is at rest) ; a comprehensive set of relations between the 
different frames can be found in Ref. Let us only notice here that the difference between the azimuthal angles 

of the final lepton in the frame of Fig. 3 and in the laboratory frame is of the 0{k\/Q'^). 

This leads to the possibility of a non vanishing transverse single spin asymmetry, the analyzing power 



which is given, according to Eqs. H51|) and (|49|l by 



^ " d6aT + dVi ' ^^^^ 



^ ^ ? da^o-^'^ z 

2 5] J d'k^ f,{x,k^) J - D'^{z,p^) 



The above equation gives the transverse SSA originated by the Sivers mechanism alone, properly taking into account 
all intrinsic motions. As it is written it depends on Xb,Q^, z^^p^ and (/)s' of course, in order to increase statistics and 
according to the experimental setups, both the numerator and denominator of Eq. (|54|l can be integrated over some 
of the variables. 

HERMES first data were gathered in the scattering of unpolarizcd Icptons (U) off "longitudinally" (L) polarized 
proton, where "longitudinal" means antiparallel to the lepton direction in the proton rest frame. Such a direction has 
a small transverse component, in the "f*p frame, with respect to the proton direction, that is 

= - smOj (1,0,0) Sj.-{Pxk±) = —sm6jSh-np sin 9^ ~ '^'m-pX^ _ y ^gg^ 

HERMES data are presented for the sin (p^ moment of the analyzing power, 

.sin 0. _ 5 Jdcj)h[da^ ~da^] sin^bh , . 

^ /d0„[daT+dai] ' ^^^^ 

which, from Eq. 1)54(1 . is 



d(f>h d'^kj_ A'^ fq/pj (x, k±) sinO^f siliip J — D'^{z,p±) sin ( 

d<P,,d'k^ fq{x,k^) J - Dl{z,p^) 



= — 7 7:rf;;^. (57) 
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where both numerator and denominator can be integrated over some of the variables. Eq. H57() gives the contribution 
of the Sivers function to A'^'^j^'^'^ . Notice that it is kinematically suppressed by the sin6'-y value, Eq. (|55|) . and that 
other contributions might be equally important; they can originate from the Collins mechanism or from higher- twist 
terms. 

More recently, data were obtained with a transversely polarized (T) proton target Sr^- {P x kj_) = sin((£— (f>s), 
and presented for the sin((/)/j — (ps) moment of the analyzing power, which singles out the Sivers contribution |lfi| . Eq. 
(|54|l in this case gives 



k, d^k± ^^fq/pi {x, k±) sm{ip - (f>s) ^ J — Dg{z,p±) sin(0ft - 



We use Eq. H58() to compute which can only receive contributions from the Sivers mechanism, and 

compare it with data, in order to gather information on the Sivers function A^/^/^t {x, k±). 



A. Parameterization of the Sivers function 

We parameterize, for each light quark flavour q = u,d, the Sivers function in the following factorizcd form p^ l32l| : 

^""U/p-^ix, fcx) = 2AA,(x) h[k^) fq/p{x, fcx) , (59) 

where 

M,{x) = N,x^^{l - xp i^^±M, , (60) 

Mfcx) = , (61) 

where iVq, Og, bq and M (GeV/c) are parameters. fq/p{x,k±) is the unpolarized distribution function defined in 
Eq. (|34|l . Since h{k±) < 1 and since we allow the constant parameter Nq to vary only inside the range [—1, 1] so that 
|A/'5(x)| < 1 for any x, the positivity bound for the Sivers function is automatically fulfilled: 

|A-/,/,r(.,..)|^^^ (62) 



'^fq/p{x,k±) 

As an alternative parameterization for h{kj_) we have also used 



h\k^) = V2e^e-'^'-/'''\ (63) 

The two parametcrizations H61|l and H63|l are indeed equivalent at low k±, but they differ at large values of kj_. 
Nevertheless, we have checked that, once multiplied by the Gaussian function of k± contained in the definition of 
f{x,k±) [see Eq. (|34|) ]. they give basically the same k± dependence to the Sivers function over the whole range, as 
shown in Fig. CKM^ = 0.25 (GeV/c)^ and M'^ = 0.36 (GeV/c)^). 
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FIG. 7: The product h{k±) ■ g{k±), solid line, compared with the product h'{k±) ■ g{k±), dashed line, as a function of fcj 



The parameterization of Eq. H63(l allows to easily perform, at 0{k±/Q), an analytical integration of Eq. H58|l . 
leading to an explicit approximate expression for the single spin asymmetry: 



.sin{tph-cj>s) 



{xB,Zh, Pt) 



(64) 



where 



Aasiv{xB,y,Zh,PT) 



2na^ 

Xb y^s 



(1 + (1 - yf) J2 4 ^^^i^-) fM D'^(^h) 



zhPt- 



/2i(fc2) 



M'(P2) (fc2 



■ exp 



ft 



(65) 



and 



where 



ao(x«, y, zn, Pt) = 2^^^ (l + (1 - vf) J2 4 f^^"^^^ ^9 (^'') 



Xb y'^s 



n{P}] 



■ exp 



It 



{kl) = 



(66) 



(67) 



and (Pt) is given in Eq. 1(211). Eq. ^ shows that A 



sin{<f>h-4>s) 
UT 



~ when Zh = or Pt = 0. 



B. HERMES data and Sivers functions 



Let us now try to understand the HERMES data on ^'J'^^*' "^^^ [lol|. according to Eq. (|58|l (exact kinematics) 
or Eqs. II64I)-(|67|I (kinematics up to 0{k±/Q)). As in Section Hill the unpolarized distribution functions are taken 
from Ref. j29l| and the fragmentation functions from Ref . |3n| | . In the numerator we take into account only the Sivers 
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N-a^ = 0.42 ±0.18 
au„ = 0.0 ±3.3 
6u„ = 2.6 ± 1.8 


iVd„ = -1.0 ± 1.8 
ffld„ = 1.1 ± 1.2 
6d„ = 5.0 ±3.6 


Nu= 1.0 ±1.9 
au = 0.52 ± 0.43 
bu^ 0.0 ±3.1 


iVj = -1.0 ±1.9 
aj= 0.0 ±4.5 
6j = 0.0 ± 2.8 


= 0.36 ± 0.43 (GeV/c)^ 


X^d.o.f. = 0.89 



TABLE I: Best values of the parameters of the Sivers functions. 



contribution of u and d quarks and antiquarks, with separate valence and sea functions. More precisely, we adopt the 
following form for the Sivers functions: 

A^/,/pT (x, k^) = 2Nq{x) h'{k^) f^/pix, k^) , (68) 
where Afq is given in Eq. (|6()|l . h' in Eq. (|63(l and q = Uy,dv,Us,ds,u,d. For the sea quark contributions we assume: 

A^/,^/pt(x,A:^) = A^/,-/pT(a;,A:^) , (69) 

for a total of 4 unknown functions, each depending on 3 parameters; in addition, h'{kj_) depends on the parameter 
M'. 

We fit the HERMES data on A™^^"^" exploiting the simplified expressions l|64|) - (|67|l . The resulting best values 
for the 13 free parameters are shown in Table 1. The errors are generated by the MINUIT minimizer. The large 
errors reflect the large errors of the data and the scarce available information. 

Our fit is shown in Fig. |H1 The solid bold line takes into account terms up to 0{k±/Q), the dashed line is 
obtained with the full exact k± kinematics, Eq. H58|l . In both cases the parameters of Table 1 are used. The 
shadowed region corresponds to one-sigma deviation at 90% CL and was calculated using the errors (Table P) and the 
parameter correlation matrix generated by MINUIT, minimizing and maximizing the function under consideration, 
in a 13-dimensional parameter space hyper-volume corresponding to one-sigma deviation. Notice that, as expected, 
the results obtained with exact or approximate kinematics are very similar. 

We show the weighted SSA A^'^'^"' as a function of one variable at a time, either Zh or Xb or Pt; the integration 
over the other variables has been performed consistently with the cuts of the HERMES experiment, at piab = 27.57 
GeV/c: 

Q2 > 1 {GcV/cf > 10 GeV^ Pt > 0.05 GeV/c 

0.023 <xs < 0.4 0.2 <Zh< 0.7 0.1 < y < 0.85 (70) 
2<Eh<l5 GcV. 

A few comments are necessary for the interpretation of the results. 

• Fig. |S1 shows that a good agreement with experiments can be obtained. However, due to the present quality 
of this first set of data, the extracted Sivers functions are not well constrained and large uncertainties are still 
possible. 

• We notice that we have checked the compatibility of the HERMES data on A^^'^'^" with the assumption 
of no Sivers effect, A^/ = 0. The data show that the probability of a zero value for the Sivers function is less 
than 0.1%. 

• It is interesting to compare the Sivers functions obtained here, Eqs. H68|) . H60|l . H63|l and Table 1, with those 
obtained by fitting the SSA observed by the E704 Collaboration in p n X processes • As stressed in 
the Introduction the question regarding the universality of the Sivers functions is a debated and open one. The 
comparison of our results with Eqs. (46)-(48) of Ref. |3 is not straightforward: one should keep in mind that 
there is no sea contribution in Ref. [3 and that the HERMES data are sensitive to much smaller x values than 
the E704 ones (which strongly depend on large x values). Moreover, the average k± and p± values adopted 
in Ref. are somewhat higher than those adopted here, derived with simplifying assumptions from data on 
azimuthal dependences in SIDIS processes. Despite all this, there is a clear indication that the two sets of Sivers 
functions arc not incompatible. The functions of Ref. 0, if used in our Eq. H58|) or H64|) - (|67|l . still allow a 
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reasonable description of the HERMES data. On the other hand, our Sivers functions of Table 1, if used to 
describe the SSA observed by E704 experiment would overestimate the data at small xp values; this could 
be easily corrected by gluon contributions (gluon Sivers function) not considered in Ref. [3 and absent at LO 
in SIDIS. 

We do not wish, at this stage and with the limited amount of available experimental information, to further 
stress such a point; the issue of the uniqueness of the Sivers functions in SIDIS and pp s-nX processes is still 
far from being phcnomcnologically established. We simply conclude that it cannot be excluded by existing data. 

• The Sivers functions obtained here are com pat ible with those extracted very recently from an analysis of Pt 
weighted HERMES data performed in Ref. [23 ■ 

Leaving aside the question of the dependence of the Sivers functions on the different physical processes, the con- 
sistency of our results can be checked within SIDIS processes, by using our functions to give predictions for other 
measured SSA. This can be done by computing, with our sets of Sivers functions, the values of A^^'^'^"''^^^ expected 
by the COMPASS experiment at CERN, which collects data in fid — ^ih^X processes at piab = 160 GeV/c, spanning 
a different kinematical region. Some preliminary results arc already available |l6l |. We neglect nuclear corrections 
and use the isospin symmetry in order to obtain the parton distribution functions of the deuterium. According to 
COMPASS experimental setup, we use the following cuts in the numerator and denominator integration of Eq. I|58|l : 

Q2 > 1 (GeV/c)2 W'^ > 25 GeV^ Pt > 0.1 GcV/c 

Eh > 4 GeV 0.2 < z^. < 0.9 O.l < y < 0.9 . (71) 

The predictions of our model are presented in Fig. 1101 and compared with the available preliminary data |l6j |. Within 
the large errors, we find a good agreement, showing the consistency of the model. 

We have also computed A^^'^'^~'^^'' for kaon production, h = K, which could be measured by HERMES; we have 
imposed the kinematical cuts of Eq. H70|l , using the fragmentation functions given in Ref. [30| . Our results are given 
in Fig. El 

Finally, we have considered the HERMES data on A^^'^'^" obtained in the semi-inclusive electro-production of 
pions on a longitudinally polarized hydrogen target 0. We have computed the Sivers contribution to this quantity, 
according to Eq. (|57|l . again with our set of Sivers functions, and compared with data. Notice that no agreement 
should be necessarily expected, as A^'^'^'" can be originated also (even dominantly) from the Collins mechanisms or 
higher-twist terms. Using the following experimental cuts: 

1 < g2 < 15 (GcV/c)2 W'^ >4: GcV^ Pt > 0.05 GeV/c 4.5 < E,, < 13.5 GeV 

0.023 <Xs < 0.4 0.2 < Zh< 0.7 0.2 < y < 0.85 , (72) 

we obtain the results depicted in Fig. 

One can conclude that our Sivers functions extracted from the HERMES data on A^^^"^"""^"^' only give a negligible 

contribution to A^'J^'^" . Not only, but the Sivers mechanism contributes with opposite signs to the transverse and 
longitudinal SSA, as can be seen from Eqs. (|57|l and H58|) . This implies that Collins mechanism and/or higher-twist 
contributions are likely to be wholly responsible for the observed , as suggested by some authors . 

V. COMMENTS AND CONCLUSIONS 

We have studied inclusive and semi-inclusive DIS processes at leading order in the QCD parton model, in the 7*p 
cm. frame and in the small Pt — Aqcd — kj_ region, where intrinsic momenta dominate the final hadron azimuthal 
and Pt distributions. We have adopted a factorized parton model scheme and exactly taken into account all intrinsic 
motions, of quarks inside the proton (fc_L) and of the final hadron with respect to the fragmenting quark {pj_). 

Wc have attempted a consistent treatment, assuming simple Gaussian k± and p± distributions and extracting from 
various sets of SIDIS data estimates about the average values {k±) and Such values are assumed to be constant, 

respectively in x and z, and to be energy independent. Simple parameterizations for the quark Sivers functions have 
been introduced. 

The resulting picture has been applied to the computation of the weighted SSA A^'^'"^" "^^^ at LO in QCD parton 
model, which directly depends on the intrinsic motions and the Sivers functions. The HERMES data clearly show 
a non zero Sivers effect; by a comparison with these data some rough estimates of the Sivers functions for u and d 
(both valence and sea) quarks have been obtained. These functions not only describe well the HERMES data, but are 
also in agreement with some COMPASS preliminary data on A^'^'^" which refer to different kinematical regions. 
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The same functions are found to give negligible contributions, with the wrong sign, to the measured longitudinal 
SSA A^"^"^" . This asymmetry can indeed be originated by the Collins mechanism and higher- twist contributions. 
Predictions for A!^^'^^~'^^'' for kaon production at HERMES have been given. 

The quark Sivcrs functions extracted from the HERMES data on pion j^^'^"' have been compared with the 
Sivers functions obtained by fitting the E704 data on SSA \x\ p it X processes. Such a comparison cannot be 
considered as conclusive, as it refers to situations with different kinematical regions and different assumptions about 
the sea contribution; however, it does not exclude the possibility that the two sets of Sivers functions - those active in 
SIDIS and in pp processes - are the same. In particular the signs seem to be the same in the two cases. Theoretical 
arguments support an opposite sign for the Sivers functions in SIDIS and Drell-Yan processes, with no conclusions 
concerning pp interactions. Our Sivers functions are compatible with those obtained in Rcf. |33 |. 

A phenomenological study of SSA and azimuthal dependences, within a factorization scheme with unintegrated 
parton distribution and fragmentation functions, is now possible. SIDIS processes with measurements of the Cahn 
effect, and the various SSA A^^jf*^ , A^^^'^'^ and A^^'^'^'^'^^^ provide a rich ground to be further explored, both 
theoretically and experimentally. 
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APPENDIX A 



It is known from symmetry principles that, within the one photon exchange approximation, the double inclusive 
cross section for unpolarizcd SIDIS processes, £p — > IhX, can have a dependence on the azimuthal angle 4'h of the 
final hadron (in the reference frame of Fig. 3) of the form 



■ ihX 



dxB dQ'^ dzh d^Pr 



^ + i? cos ( 



Ccos2(j)h 



(Al) 



where A, B and C are scalar quantities which do not depend on 0^ . This is explicitly visible in the approximate 
expression (|38|l and we wonder whether Eq. (|31|) satisfies in general such a condition. 
From Eqs. 1(23, (E!)) and (|2H1) one can see that Eq. (EU is of the form 



J d^kx_ [a + bt-k^ + c{t-kx_f] F{Pt • . . .) 



(A2) 



where a, b and c do not depend on angles and the . . . stands for scalar variables which also do not depend on azimuthal 
angles. 

As a consequence, a tcnsorial analysis of Eq. ljA2|) shows that Eq. H31|) can only contain azimuthal dependences 
through the integrals: 



cos <ph 



be- J d^k±k±F{PT -k^) ^ e- F 

j d^k^ ik^)j F{Pt ■ fc±) ~ ^^ (Prh {Pt)j ~ cos^ 0„ = ^ + 



C ^1 



which agree with Eq. (|A1|I . 
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FIG. 8: HERMES data on A 
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Hol for scattering off a transversely polarized proton target and pion production. The 
curves are the results of our fits, with exact kinematics (dashed line) or keeping only terms up to 0{k±/Q) (solid bold line). 
The shadowed region corresponds to the theoretical uncertainty due to the parameter errors. Data on the Pt dependence of 
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are not available yet, and the third column of the figure gives our predictions. 
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FIG. 9: Predictions for kaon asymmetries A^^'*''"'*^' at HERMES kinematics for scattering off a transversely polarized proton 
target. The curves correspond to calculation with exact kinematics (dashed line) or keeping only terms up to 0{k±/Q) (solid 
bold line). The shadowed region corresponds to the theoretical uncertainty due to the parameter errors. 
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FIG. 10: COMPASS preliminary data ^ on A'^^^!'^''"'^^' 

production of positively (/i'*') and negatively (h~) charged hadrons. The curves show our predictions, according to the values 
of the parameters for the Sivers functions given in Table 1 and obtained from fitting the HERMES data on J^^^r^"^^'"^^^ . Again, 
the dashed line refers to exact kinematics, Eq. 1581 . while the solid bold line is obtained by keeping only terms up to 0(k^/Q). 
The shadowed region shows the theoretical uncertainty due to the parameter errors. 



21 



^ 0.075 
< 0.05 



0.025 


-0.025 
^ 0.075 

< 0.05 
0.025 



-0.025 
^ 0.075 

< 0.05 
0.025 



-0.025 



71 



71 



71 



71 



71 



71 




71 




0.5 1 

Pj (GeV/c) 



FIG. 11: HERMES data on A'^'^'^^ for scattering off a longitudinally polarized proton target and pion production. The 
curves show the contribution of our Sivers functions alone, with exact kinematics (Eq. 1571 , dashed line) or keeping only terms 
up to 0{k±/Q) (solid bold line). The shadowed region corresponds to the theoretical uncertainty due to the parameter errors. 



